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Abstract. Recent simulations of supernova-driven turbulence within the ISM support the ex- 
istence of a large-scale dynamo. With a growth time of about two hundred million years, the 
dynamo is quite fast - in contradiction to many assertions in the literature. We here present 
details on the scaling of the dynamo effect within the simulations and discuss global mean- 
field models based on the adopted turbulence coefficients. The results are compared to global 
simulations of the magneto-rotational instability. 
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1. Introduction 



Lar ge-scale coh erent magnetic fields are observed in many spiral galaxies (jBeck et al 



1996; Beck 2009). It is now widely accepted that a turbulent dynamo process is re- 



sponsible for the sustained amplification of the mean magnetic field within galaxies. The 
turbulence driven by supernova explosions together with the diffe rential rotation should 



in fact lead to the action of the classical afl dynamo (jParker 11971 ). First s teps to globall y 



model flat objects in an ellipsoidal geometry were done bv IStix 



( 197,^ : IWhiteT(|l978h . 



Starting from the early nineties, fir st multidimensional numerical mean- field models as - 
suming a disk geometry appeared ( Donner fe Brandenburg I Il990l ; Elstner et al. 1992 ). 
These models applied simple estimates for the turbulent transport coefficients in the 
m ean-field induction equat i on. 

Riidiger fc Kitchatinov (1993), in the framework of second order correlation approx- 



imation (SOCA), derived th e full a tensor and turbulent diffusivity for stratified tur- 
bulence. Froh lich fc Sch ultz I (|l996l) . under the assumption of hydrostatic equilibrium, 
computed the turbulent velocity distribution u'(z) based on the observed density strat- 
ification and an observationally motivated gravitational potential from stars and dark 
matter. Applying SOCA theory, one could now specify the full a tensor, resulting in a 

model with the correlation time t c of the turbulence as the only free param eter. 

With the increased computer power of the nineties. IKaisig et al. ( 19931 ) and IZiegler 



1996), in an alternative approach, computed the turbulence coefficients by averaging 



the results from simulations of a single supernova remnant. It, however, turned out that 
the values from isolated explosions were about one magnitude smaller then the estimated 
numbers from the pre vious considera tions. Similar results, based on semi-analytical mod- 
els, were obtained bv iFerrieTel 1 19981 ). Extending the investigation to super-bubbles and 
including the stratification of the ISM, she finally found turbulence coefficients of ex- 
pected magnitu de. However, with a huge es cape velocity, which could switch off the 
dynamo process. iFerriere and Schmitt 1 (|2000h used these results as input for an axisym- 
metric dynamo model of our Galaxy. They reached growth times of about 450 Myr but 
with a ten times larger toroidal field than the radial component, i.e., with a rather low 
pitch angle. 

A lot of models, but mostly axisymmetric, were published during that period. Most 
of the models were investigating different aspects of the turbulence, looking also for 
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parameters where non-axisymmetric solutions could appear. Another application was the 
effect of radial flows on the dynamo properties driven by the Maxwell stress of the dynamo 
generated magnetic field. Models inclu ding a wind were investigat e d, especially to explai n 
vertical fields above the disk - notably ferandenburg et all ( 19931 ); Elstner et al. I ( 19951 ). 
Both models use d different assumptions abo ut the halo diffusivity and the wind structure. 
In the model of Branden burg et al. I (|l993l) . diamagnetism was included, and because 
of the positive turbulence gradient, the pump ing was not an escape velocity but an 
inward transport term. In Elst ner et al. I (|l995l ). a passive halo with low conductivity 
was assumed. These models already showed a fast growth time below 0.1 Gyr. All these 
models could in principle explain the main features of the observed magnetic fields in 
nearby spiral galaxies. But the connection between optical and magnetic spirals was still 
unclear. 

With the advent of the new century, full-blown 3D models became feasible, and one 
could now attack the spiral structure. Simple artificial velocity models adopted from 
density- wave theory were now taken into account. One specta cular observational resul t 
was the anti-correlation of optical and magnetic arms, found bv lBeck &: Hoerries" ( 1996t ). 
A turbulent dynamo could reproduce such a behaviour due to the diff erent properties of 
the tu r bulence in t he sp iral arms and in the int er-arm regions - see 
1 1998MShukurov I ljl998f) and lRohde etaD (Il999h . 

More elaborated models used the v elocity from numeric al galaxy models computed 
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with SPH or sticky-particle methods ( Elstner et al. 2000). Models for spiral galaxies 



produced no inter-arm fields by the large-scale spiral motion alone. Only with a variation 
of the t urbulence, one could have inter-a rm fields. The situation was different for bar 
galaxies ( Otmiano wska-Ma zur et al. 2002). In these simulations, one could observe inter- 
arm magnetic fields in the outer spiral ar m, and the field la gged behind the optica l 
arm (jMoss et al. I Il998l cf.). More recently, iBeck et al.l (|2005l ) and iMoss et~aTl (|2007l ) 
conducted models for barred galaxies with stationary flows. 

Another problem was raised in the context of the so-called catastrophic quenching: 
Because of helicity conservation, in the ideal MHD case, the a term will be quenched (with 
the square-root of the magnetic Reynolds number) if there is no additional transport of 
magnetic helicity in the box. Under the assumption of a galactic wind or fountain flow, 
these proble ms may disappea r. But a too strong wind would also result in a suppression of 
the dynamo (|Sur et al.ll2007l ). Recent investigations have shown that already the galactic 
shear flow can supp ress the catastrophic quenching of the dynamo (jBrandenburg I [20091 : 
Kapyla et al~ll2008h . 



The rece nt box model of th e turbulent ISM driven by multiple clustered supernova 
explosions (Gresse l et al. 2008), for the first time demonstrated the presence of a dynamo 
process by direct numerical simulations. It turned out that the turbulence coefficients 
derived from the SOCA models are in good agreement wit h the new resu l ts fro m the 



direct simulations. In a similar way, the box simulations of lHanasz et al. I (|2004f ) have 
shown an exponential growth of the magnetic field. Their simulations, however, do not 
resolve the supernova remnant in itself. Instead, the explosion site only serves as a source 
for the released cosmic ray energy. Due to the fast diffusion of cosmic rays, the pressure 
redistribution is fast enough to allow for similarly short growth times of several 100 Myr. 
The role of a magnetic instability of the Parker type could be essential for this process, 
too. 

In general, the effect of magnetic instabilities in releasing kinetic energy should be 
considered as a possible mechanism for magnetic field amplification in galaxies - espe- 
cially in regions with low star formation activity. The primary candidate of this type 
is the magneto-rotational instability (MRI), which has a growth time on the order of 
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the galactic rotation frequency and thus can become efficient on very short timescales. 
With todays computational resources, global mode ls of the MM are already within reach 



( Dziourkevitch et al. 2004 ; Nishikori et al. 2006f h However, these models still have to 



assume quite strong initial fields to properly resolve the fastest growing modes of the 
instability. 



2. The pitch angle problem 

When it comes to the explanation of the observed magnetic fields in spiral galaxies, 
there is one main reason for favouring the dynamo hypothesis over the primordial field 
hypothesis, namely the existence of fields - of equipartition strength - with a large 
pitch angle. If the field amplification were mainly due to the shearing motion of the 
differential rotation, the field would have to be wound up several times to explain its 
observed strength. This, however, would lead to an amplification of the toroidal field 
only, resulting in a configuration with an intrinsically small pitch angle. 

For the dynamo scenario, which furthermore implies exponential growth, we addition- 
ally have an amplification of the poloidal field. The pitch angle p of the growing field in 
the kinematic regime of the aQ dynamo can be estimated from the dimensionless dynamo 
numbers C a and Co (for a definition see below) via 

(2.1) 

This leads to values between 5° and 30°, for reasonable values of C a = 5-10 and Cn = 25 
800. For a given rotational velocity, the angular frequency is defined by the turnover 
radius, i.e., the point where the rotation curve becomes flat. Small turnover radii lead 
to rather large Cn and thus to low pitch angles. In the saturated regime, the situation 
becomes worse - at least for the conventional aQ dynamo. This is because a will then 
be reduced by non-linear effects and the pitch angle fades. 

Direct simulations of the ISM turbulence show a possible way out of this dilemma, 
namely by another saturation process, which we will demonstrate in the following sec- 
tion. Moreover, the action of the magneto-rotational instability may be a reasonable 
candidate for the explanation of large pitch angles as they are observe d in maps of po 




larised synchrotron emission for instance in the ringed galaxy NGC 4736 (IChvzv fc Buta 
20081) . 



3. Global models with turbulence coefficients from direct simulations 

Global simulations of interstellar turbulence remain out of reach, at least if one is 
interested in resolving the relevant scales. Therefore, we apply a hybrid approach, where 
we consider the mean induction equation 

OB -* — -* 

= curl(w xB + aoB — t]to curl_B) (3-1) 

for the mean magnetic field B and apply closure parameters obtained from direct simu- 
lations. The background rotation serves as the main part of the mean flow 

(3.2) 
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Figure 1. Radial dependence of Q and vertical profiles of the fountain flow and turbulent 

transport coefficients. 

and is supplemented by a vertical wind u z = uqz which is a result of the density- 
stratified turbulence driven by SNe (cf. Fig. 2 in Gressel et al. 2009). The test-field 
method, furthermore, yields an a tensor of the form 



(3.3) 



where the dependence of these quantities on the height z is modelled according to the 
profiles in lGressel et al. I ( 20081 ). 

For a realistic global galaxy model one further needs the radial dependence of the 
relevant effects. The radial profiles of the rotational frequency, SN rate, gravitational 
potential and density are in principle known. By means of a parameter study, we try 
to obtain the a effect from our box simulations, accordingly. From these, we find a 
dependence of the diagonal terms of the a tensor on rotation, namely cv w oc fl 5 . 
This is weaker than expected from SOCA theory, which predicts a scaling linear with 
Q and only a minor rotational quenching. Consistent with SOCA, we find the turbulent 
diamagnetism ^y z to be nearly independent of f2. For the turbulent diffusivity r) t , we apply 
only horizontal components which also vary according to f2 0,5 . The dependence on the 
supernova rate is given in Gressel et al. (2009). 

After proper normalisation with a length of I = 1 kpc and a diffusion time = I 2 /rjo 
with 770 = 1 kpc 2 Gyr" 1 , the dynamo is characterised by the Reynolds numbers 



C a = a l/r/o , C n = fl l 2 /i] and C w = u l/r) . 



(3.4) 



We assume a disk of radius R = 5 kpc and with a height of H = 2 kpc, according to our 
box dimension in the direct simulations. We approximate the profiles of the turbulence 
coefficients by smooth functions (see Fig.[T]). The vertical dependence of the a tensor can 
simply be reproduced by a sine function 



a(z) oc sin(z7r/h Q ) 



(3.5) 



with a scale height h a — 1.5 kpc. The mean vertical velocity is linearly growing with 
height. Only the diagonal terms of the a tensor have a radial dependence proportional 
to whereas the other terms are constant. We, furthermore, use a simple a quenching 



ao/(l + B 2 /B 2 eq ) 



(3.6) 



for the full tensor, i.e., including the diamagnetic term j z . 

In our first run, we neglect turbulent diamagnetism and mean vertical velocity, as it is 
done in the classical afl dynamo. In this model, we find a growth time of 700 Myr, a field 
strength for the toroidal component B v = 20B cq and the radial component B r — 0.6B cq . 
If we neglect only the fountain flow, the field decays. For the real model, we find a 
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Figure 2. Left: Time evolution of magnetic energy for model A, with wind (dashed) and model 
F, without wind and pumping (solid). Middle: Time evolution from 1.9 Gyr until 9.6 Gyr of 
the radial pitch angle distribution for model F. Right: Pitch angle for model A at 0.5 Gyr and 
9.1 Gyr. 
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Figure 3. Magnetic vectors as observed in polarised emission (neglecting Faraday rotation) 
for the final fields in model F (left) and model A (right) 

growth rate of 250 Myr - similar to the box simulation. The field now saturates earlier 
with B v = l.lB cq and B r = 0.3-B cq . 

The final state is not determined by the usual a quenching in this case. Instead, 
a combination of fountain flow and quenching of the diamagnetic term leads to the 
saturation of the dynamo. This means, the dynamo stops growing because of field losses 
by the fountain flow. As a nice consequence of this non-linear feedback, the pitch angle 
remains large in the saturated regime (see Fig. 

Moreover, in Fig. [3]we present directions of the magnetic field resulting from polarised 
emission of synchrotron radiation, where we assumed a scale height of 500 pc for the 
relativistic electron distribution. 

In Table [TJ the pitch angle is given for a choice of different models. Comparing p = 
arctan( 1k /C ct /Cn) with the initial pitch angle po during the kinematic growth phase, we 
find reasonable agreement. The final value pfi na i in the saturated state for the first three 
models is similar to the value during the kinematic growth phase. 

In order to have a higher pitch angle for our rotation law with tq = 1 kpc, we need 
larger values of C a and, consequently, a larger vertical outflow. Otherwise the a quenching 
would reduce the final pitch angle - as is the case in model E. This increase of the mean 
outflow velocity with increasing SN rate is indeed observed in the box simulations (cf. 
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Table 1. Pitch angles for different dynamo models 
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Gressel et al. 2009). For model F, one can see a dramatic reduction of the pitch angle 
due to the missing mean vertical velocity and diamagnetic pumping. 



4. Magneto-Rotational Instability 

Galaxies are known to be u nstable with resp ect t o the magneto-rotationa l instability. 
This instability, discovered by IVelikhov I (jl959h and Chandrasckhar I (Il960h. has mean- 
while been studied in great detail in t he astrophysical cont e xt (se e e.g. Balbus fc Hawlevl 
199lt iKitchatinov fe Mazur I Il997t ). ISellwood fc Balbus I (|l999h pointed to the possi- 



bility of turbulence genera tion by the MRI in regions of low star formation activity. 



e get 

Piontek fc Ostriker I (|2007f ) recently simulated the MRI in a vertically stratified two- 



p hase ISM box model. 

Kitchatinov fc Riidiger (2004) investigated the stability of global galactic disks. For 
a disk threaded by a weak vertical field, they estimated a field strength of 10~ 25 G for 
the onset of the instability. Accordingly, the MRI should already be present during the 
galaxy formation process - and therefore a good candidate for the seed field of a turbulent 
dynamo. The most unstable mode in their idealised galactic disk model is of quadrupolar 
type. A very interesting result is the large pitch angle of 45° for the fastest growing mode. 
For weak fields, the unstable modes have a small spatial scale but with the property of 
a uniform large pitch angle. This leads to a smooth map of magnetic vectors for the 
polarised radio emission, but the Faraday rotation measure will be small for these types 
of magnetic fields. 

The instability, moreover, has a very short growth time proportional to the rota- 
tion time. A proof for a large scale dynamo driven by the MRI in the sense, that a 
small initial field can be amplified over several magnitudes has not been given yet. This 
is because it is numerical difficult to resolve the small spatial scales for the unstable 
modes connected to the weak initial field. Nevertheless, there are hints that the process 
is possible. Especially the turbulence excited by MRI in a stratified disk gives a non- 



vanishing a effect, but with negative sign in the northern hemisphere (jArlt fc Riidiger 
200lUBrandenburg fc Sokoloff Il2002h . 



Since we are interested in configurations with large pitch angles, we here consider 
global models of a density-stratified galactic disk with a uniform vertical field. Tests 
with toroidal fields did not produce a comparable radial field component. For simplicity, 
we assumed an isothermal disk where the rotation and vertical density stratification is 
supported by the gravitational potential. Without magnetic fields the configuration was 
tested to be in a stationary equilibrium over 10 Gyr with only small velocity fluctuations 
below 0.01 kms -1 , which are slowly decaying with time. 

The applied density profile is 
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with the scale height h p = 0.3 kpc and pa — 10 24 gcm 3 . Note that the density minimum 
is about 0.02 pa. The rotation curve is, furthermore, defined by (|3.2p with the turnover 
radius r n = 0.35 kpc and fl Q = 200 Gyr" 1 . The inner radius of our disk model is r ln — 
lkpc and the outer radius r out = 4 kpc. The vertical extent is ±5 kpc. Within ±lkpc, 
we use a constant grid spacing with 256 grid points and a non-uniform grid for the outer 
regions with 128 points at each side. This means, the inner disk has a resolution of about 
7 pc and for the outer disk we yield about 80 pc. This is consistent with a density variation 
of a factor 100, in order to get a comparable growth of the unstable modes for a uniform 
magnetic field strength within our restricted resolution. 

The radial boundary is reflecting and the vertical boundary is reflecting for the velocity 
only, while the tangential component of the magnetic field and the gradient of the vertical 
component are set to zero. Under these conditions, the vertical flux is conserved but a 
development of horizontal flux is allowed. Applying a weak and uniform vertical magnetic 



V. at time=1.78 



V. at time=4.03 




Figure 4. Radial profiles of tp, z averaged velocity components in [kms -1 ] (first row), rms-ve- 
locity components in [kms -1 ] (second row), ip, z-averaged magnetic field components in [pG] 
(third row) and rms magnetic field components in [pG] (fourth row) at 1.8 Gyr 



field of 0.1 pG, we find an exponential growth of the unstable mode after a transient 
relaxation phase of about 200 Myr. The field starts to grow at the inner radial boundary, 
where fl has a maximum. As the field grows further, it spreads out radially. At the same 
time, the Maxwell stress causes an inflow in the inner part of the magnetised region. 
At the outer magnetic region, we observe an outflow of mass, which is a consequence of 
angular momentum conservation. 

Because of our closed boundaries, the density increases strongly at the inner boundary 
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and moderately at the outer boundary. We see a weak change of the rotation curve at the 
boundaries, compensating the developing pressure gradients. This is probably an artefact 
of the closed boundaries. Tests with larger radial domains have shown that the change 
of the rotation curve is restricted to the boundary region. 

Finally, the field growth saturates. At the moment it is unclear whether the saturation 
is caused by the non-linear feedback of the developed turbulence. The reduction of the 
pitch angle could in fact be a hint that the system reaches the diffusive limit of the 
instability, where we know from linear analysis that the unstable mode, in this regime, 
is mainly toroidal. On the other hand, we see a large-scale meridional circulation in the 
disk. Moreover, our model disk has a limited mass reservoir, which could be another 
reason for the early saturation with rather low magnetic field strength. As long as the 
toroidal and radial field strength are below the field strength of the initial vertical field, 
the field parallel to the disk grows with a time of 0.05 Gyr; later we observe growth 
times of 0.8 Gyr until the saturation sets in. It should be noted that the growth-time of a 
possible MRI-driven dynamo has not to be as fast as the growth time of the instability. 

In Fig. 21 we plot radial profiles which have been averaged along the vertical and 
azimuthal direction for —5 kpc < z < kpc and < ip < 2ir. The values shown are the 
components of the velocity, magnetic field and their corresponding rms values. A large- 
scale meridional circulation develops with an amplitude of roughly 0.5 of the rms value 
(~ 6kms~ 1 for the vertical component and ~ 3kms _1 for the radial component). 

The mean magnetic field along the z and tp direction are of the same order as the 
rms values of the fields. Only the mean radial field is by an order of magnitude smaller 
than the rms field. We find B r r ms / ' B™ s = 0.3 and for the mean-field B r /B v = 0.05 - 0.1, 
corresponding to pitch angles of 17° and 3° — 6° , respectively. In contrast, the pitch 
angle in the polarisation map varies from 10° to 20° if the boundaries are excluded 
(1.3 kpc < r < 3.6 kpc). 

Initially, at time 1.8 Gyr, the pitch angle varies from 20° to 40° for 1.3 kpc < r < 2 kpc. 
In the outer region the field is still too weak, and we do not observe any polarised emission. 
Fig. [5] shows a slice through the magnetic field at z = —0.8 kpc, the vertically averaged 
field along the lower disk half (—5 kpc < z < Okpc) and the magnetic direction of 
polarised emission (neglecting Faraday rotation). Varying the scale height height of the 
relativistic electron distribution from 0.5 kpc to 2.0 kpc gave no significant differences. 
Note the smooth polarisation map for an incoherent magnetic field (i.e., with many 
reversals along the vertical direction) . The pitch angle is here much larger than the pitch 
angle for the averaged field. 

5. Summary 

New direct simulations of supernova-driven turbulence in a Cartesian box have shown 
the existence of the turbulent dynamo in a galactic environment. We find the resulting 
balance of the turbulent diamagnetism with the mean vertical outflow to be crucial for 
the operation of the dynamo. Using the turbulence of the direct simulations for the 
determination of the transport coefficients, one can show that the fast dynamo growth 
time of a few hundred million years is due to the fast transport by the mean vertical 
fountain flow and the downward pumping of the turbulent diamagnetism. This is in 
contrast to the slow conventional att dynamo with growth times of the order of 1 Gyr, 
where only diffusion contributes to a vertical field transport. 

We have shown the possibility of an alternative saturation process other than the 
simple quenching of the diagonal a term. In the case of quenching via the transport 
mechanisms, the pitch angle of the final saturated field will not deplete. But still, there 




Figure 5. Magnetic field slices at z — — 0.8kpc (first row) and averaged field over one disk half 
(second row) with colour coded z component. Magnetic field direction from polarised emission 
with colour coded total intensity (third row). 



are cases of very strong differential rotation where it may be difficult to get observed 
pitch angles of up to 40° with actual SN-driven dynamo models. 

More detailed investigations of the ISM model must be performed. For instance, cos- 
mic rays could further improve the situation. Also the effects of strong magnetic fields 
have to be considered, as magnetic instabilities could contribute to further amplification 
processes. The magneto-rotational instability would be a good candidate - not only as a 
generator for the seed field, but also for the generation of magnetic fields with large pitch 
angles, ft is not clear, however, whether an MR1 driven dynamo can generate enough 
vertical flux through the disk in order to get solutions with large enough pitch angles, 
which are known from the linear unstable modes. But at least under a given vertical flux, 
the MR1 leads to incoherent magnetic fields, which appear in the polarised synchrotron 
emission (neglecting Faraday rotation) as smooth spiral fields. The mean-field generated 
in such models has a smaller pitch angle than the polarised emission vectors, so a signif- 
icant contribution from a small scale field enhances the pitch angle in the polarisation 
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map. Future RM synthesis observations would be very helpful for the understanding of 
magnetic fields in external galaxies. 
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